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LITERATURE  REVIEW  OF  MARINE  WETLAND  AND  ESTUARINE 
WATER  QUALITY  AMD  ECOSYSTEM  MODELS 

PART  I :  INTRODUCTION 

1.  This  literature  review  surveys  published  models  of  hydro¬ 
dynamics,  sediment  transport,  water  quality,  and  ecosystems  in  estuaries 
and  intertidal  wetlands.  This  review  was  based  on  the  needs  and  prob¬ 
lems  of  the  U.  S.  Army  Corps  of  Engineers  (CE)  in  the  coastal  zone  as 
revealed  by  a  survey  of  CE  Field  Offices  with  coastal  responsibilities. 
In  addition,  a  workshop,  conducted  in  New  Orleans  in  June  1979 >  assess¬ 
ing  the  state  of  the  art  in  understanding  and  modeling  estuarine  and 
wetland  processes,  served  to  define  the  scope  of  the  literature  review. 
Both  the  workshop  report  and  the  report  of  the  survey  of  CE  Field 
Offices  have  been  published  (Hamilton  1979a,  1979b). 

2.  Activities  of  the  CE  primarily  alter  the  physical  system  in 
estuaries  and  wetlands  by  altering  freshwater  inflow,  dredging  naviga¬ 
tion  channels,  constructing  harbors  and  shore  protection  structures, 
and  creating  intertidal  marsh.  Therefore,  prediction  of  the  effects 
of  these  activities  on  the  physical  regime  is  of  primary  importance; 
the  hydrodynamic  models  that  are  the  basis  of  the  majority  of  predic¬ 
tive  techniques  are  surveyed  and  discussed  in  Part  II  of  this  report. 
Alteration  of  the  physical  regime  can  also  profoundly  affect  sediment 
transport,  biological  and  chemical  processes,  and  possibly  fisheries. 
Thus,  the  difficult  problem  of  modeling  sediment  transport  in  estuaries 
is  reviewed  in  Part  III.  Part  IV  concerns  the  important  problem  of 
modeling  water  quality  parameters  in  estuaries.  The  majority  of  models 
considered  are  based  on  mass  conservation  that  requires  the  use  of 
hydrodynamic  submodels  to  predict  dispersion  and  velocity  fields.  How¬ 
ever,  some  water  quality  models  based  on  population  dynamics  and/or 
energy  flow  are  also  considered  because  of  possible  application  to 
toxicity  effects  on  estuarine  organisms.  Some  of  the  more  sophisti¬ 
cated  water  quality  models  of  estuaries  are  effectively  applied 
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ecosystem  models  of  the  lower  trophic  levels.  Ecosystem  models  of 
estuaries,  which  have  usually  developed  from  marine  models,  are  reviewed 
in  Part  V.  Finally,  the  much  less  developed  field  of  marsh  ecosystem 
modeling  is  discussed  in  Part  VI. 

3.  This  review  has  been  written  as  a  guide  to  the  published  lit¬ 
erature  on  estuarine  and  wetland  models  and  attempts  to  span  the  gap 
between  a  simple  reference  list  and  a  comprehensive  critical  review. 
Important  studies  of  the  last  five  years  are  put  in  historical  perspec¬ 
tive,  and  there  is  some  discussion  of  their  major  achievements  and 
limitations. 

k.  The  authors  of  this  survey  have  attempted  to  highlight  those 
model  studies  of  estuaries  and  wetlands  that  are  or  have  been  important 
contributions  of  the  state  of  the  art.  The  authors  take  full  responsi¬ 
bility  for  their  views  and  prejudices  and  sincerely  apologize  if  impor¬ 
tant  models,  which  are  relevant  to  this  investigation,  have  been 
inadvertently  overlooked. 
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PART  II:  HYDRODYNAMIC  MODELS  OF  ESTUARIES 


Introduction 


5.  The  basis  of  all  water  quality  models  and  many  ecosystem  mod¬ 
els  that  use  the  principal  of  conservation  of  mass  is  a  model  of  the 
circulation  of  the  estuary  or  coastal  bay.  The  investigation  of  hydro¬ 
dynamics  of  estuaries  by  the  means  of  models  has  a  long  history.  Most 
of  the  earlier  studies  concerned  the  propagation  of  tidal  waves  in  chan¬ 
nels  (Defant  1961,  Dronkers  196U)  and  only  recently  has  the  internal 
depth  dependent  circulation  been  a  subject  of  model  studies. 

6.  The  circulation  in  a  partially  mixed  estuary  or  bay  is  mainly 
driven  by  the  tide,  river  inflow,  salinity  gradient,  and  wind.  The 
importance  of  wind  has  only  recently  been  realized.  Based  on  extensive 
field  observations,  particularly  long-term  current  meter  measurements  in 
Chesapeake  Bay  (Wang  1980),  it  is  now  thought  that  in  many  estuaries 
all  four  forcing  mechanisms  are  equally  important  in  influencing  the 
circulation  and  the  transport  of  pollutants  and  sediments.  Estuarine 
circulations  are  perhaps  some  of  the  most  complex  geophysical  flows  in 
existence.  Besides  the  complexity  of  the  driving  forces,  there  is  the 
difficult  phenomena  of  turbulence  in  a  stratified  system  and  the  influ¬ 
ence  of  channel  topography  in  inducing  cross-channel  circulations  and 
hence  cross-channel  variability.  A  review  of  recent  turbulence  studies 
in  highly  stratified  estuaries  is  given  by  Gardner  et  al.  (1979)  who 
show  that  significant  mixing  can  arise  from  tidal-induced  internal  waves 
or  bores  as  well  as  due  to  boundary  layer  or  wind-induced  turbulence. 

A  more  traditional  review  is  given  by  Bowden  (1977)  who  emphasizes  the 
need  for  parameterizing  the  vertical  eddy  viscosity  and  diffusion  coef¬ 
ficients  in  order  to  improve  models  of  vertical  circulation.  Topo¬ 
graphic  effects  have  had  little  study.  However,  Dyer  (1977)  discusses 
possible  secondary  circulations  induced  by  changes  in  channel  topography 
and  an  extensive  field  experiment  has  been  performed  at  a  bend  in  the 
lower  Potomac  Estuary  to  investigate  these  effects  (Boicourt  1979). 

7.  Physical  oceanography  of  estuaries  is  the  subject  of  a  number 


of  texts  by  Dyer  (1973),  Officer  (1976),  and  McDowell  and  O’Connor  (1977) 
and  also  symposia  edited  by  Kjerfve  (1978)  and  Hamilton  and  Macdonald 
(1980).  For  a  review  of  numerical  models  of  estuaries  and  coastal  bays, 
see  Liu  and  Leendertse  (1978). 


One-Dimensional  Models 


8.  Apart  from  tidal  propagation,  one  of  the  earlier  uses  of  one¬ 
dimensional  models,  which  include  an  advection-diffusion  equation  for 
a  solute,  is  the  investigation  of  salinity  intrusion  in  an  estuary 
(Thatcher  and  Harleman  1972;  Stigter  and  Siemons  1967 ;  Williams  and 
West  1973).  Neglect  of  the  variation  with  depth  of  current  and  salinity 
means  that  the  important  upstream  flux  of  salt  due  to  density  currents 
as  well  as  vertical  phase  differences  in  the  tidal  variations  of  current 
and  salinity  are  all  lumped  into  one  parameter:  the  one-dimensional 
longitudinal  dispersion  coefficient.  It  is  not  surprising,  therefore, 
that  Williams  and  West  (1973)  found  that  the  dispersion  coefficient  is 
highly  dependent  on  riverflow  and  distance  from  the  estuary  mouth. 

These  kinds  of  one-dimensional  models  are  difficult  to  use  in  a  predic¬ 
tive  mode  in  some  partially  mixed  estuaries.  These  models  are  also  the 
hydrodynamic  basis  for  the  majority  of  water  quality  models  of  estuaries. 

Two-Dimensional  Models 


Longitudinal  and  vertical  dimensions 

9.  The  inclusion  of  the  depth  dimension  makes  possible  the  inclu¬ 
sion  of  the  majority  of  the  fundamental  physical  process  that  contributes 
to  advection  and  diffusion  of  salt  and  pollutants.  The  first  study 
which  showed  the  importance  of  the  depth  variation  of  current  and  salin¬ 
ity  was  the  steady-state  analytical  solution  of  Rat'tray  and  Hansen  (1962). 
This  was  later  expanded  to  include  the  horizontal  salt  balance  explicitly 
by  Hansen  and  Rattray  (1965).  A  similar  approach  was  also  used  for 
fjords  by  Rattray  (1967)  and  Winter  (1973).  The  same  set  of  steady- 
state  equations  was  solved  numerically  by  Festa  and  Hansen  (1976). 
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10.  The  time-dependent  equations  for  a  tidally  driven  partially- 
mixed  estuary  were  solved  numerically  by  Hamilton  (1975)  and  Blumberg 
(1975)  and  applied  to  the  Rotterdam  Waterway  and  the  Potomac  Estuary, 
respectively.  Both  numerical  methods  were  explicit  (i.e.  time  step 
limit  by  the  gravitational  wave  speed — the  Courant-Friedrich-Lewys 
condition)  and  had  small  steps  of  approximately  1  min.  Semi-implicit 
extensions  which  enabled  time  steps  of  about  one  half  hour  to  be  used 
with  little  loss  in  accuracy  have  been  developed  by  Hamilton  (1976)  and 
Wang  and  Kravitz  (1980).  This  allows  the  important  subtidal,  wind-driven 
transients  to  be  investigated  efficiently  (Wang  1980).  These  kinds  of 
models  have  also  bfeen  applied  to  branching  estuaries  (Elliott  1976). 

11.  The  important  parameters  in  these  vertical  circulation  models 
are  the  vertical  eddy  coefficients.  The  effects  of  different  formula¬ 
tions  of  these  coefficients  are  shown  in  Bowden  and  Hamilton  (1975),  and 
the  importance  of  the  effects  of  stratification  on  vertical  turbulent 
diffusion  is  stressed.  The  turbulent  stresses  may  be  calculated  from 
the  conservation  equation  for  turbulent  kinetic  energy  with  seme  assump¬ 
tions,  and  this  approach  has  been  used  for  an  unstratified  tidal  channel 
by  Johns  (1978).  This  seems  to  be  a  promising  approach  to  overcoming  the 
limitations  of  first  order  K-theory  to  parameterize  the  eddy  coeffi¬ 
cients,  but  further  complications  are  introduced  such  as  very  fine 
grids  to  resolve  the  boundary  layers. 

12.  The  models  discussed  above  consider  the  depth  coordinate  to 
be  continuous  and  in  the  numerical  models  this  is  approximated  by  a 
fixed  grid.  However,  fjords  such  as  Puget  Sound  and  salt  wedge  estuaries 
such  as  the  Mississippi  are  highly  stratified  and  often  more  appropri¬ 
ately  modeled  by  using  layers.  The  displacement  of  the  interfaces  is  a 
function  of  space  and  time.  A  time-dependent,  two-layer  model  of  a 
deep  fjord  (Alberni  Inlet,  British  Columbia)  has  been  developed  by 
Hodgins  (1979)  and  for  a  salt  wedge  estuary  by  Vreugdenhil  (1970). 
Horizontal  dimension 

13.  The  vertically  integrated  equations  of  motion  have  been  used 
for  storm  surge  modeling  of  coastal  seas  and  estuarine  bays  for  a  number 
of  years.  Reviews  of  the  different  kinds  of  models  are  given  by  Hinwcod 
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and  Wallis  (1975a,  1975b).  The  use  of  these  models  to  predict  water 
level  and  depth  averaged  current  variations  due  to  tides  and  wind  is 
now  quite  sophisticated  in  that  complex  geometry  using  variable  grids 
and  imbedded  subgrid  scale  features  such  as  channels,  barriers,  and 
canopies  can  now  be  handled  (Butler  1979;  Reid  et  al.  1977a,  1977b). 

A  good  example  of  application  of  this  type  of  model  to  tidal  prediction 
in  Chesapeake  Bay  is  given  by  Blumberg  (1977);  however,  there  are  many 
examples  of  similar  applications  to  bays  in  the  literature  (Gordon  and 
Spaulding  197*0. 

1*+.  All  the  numerical  models  mentioned  so  far  in  this  section 
employ  finite  difference  methods.  Models  using  finite  element  methods 
are  based  on  the  minimization  of  functions  defined  over  elements  (often 
triangular)  subject  to  the  constraints  of  the  equations  of  motion.  For¬ 
mulation  of  these  models  is  generally  more  complex  than  their  finite 
difference  counterparts  and  less  is  known  of  their  stability  and  accuracy 
behavior.  Examples  are  studies  by  Pearson  and  Winter  (1977),  Wang  and 
Connor  (1975),  and  Gray,  Pinder,  and  Brebbia  (1977). 

15.  The  addition  of  a  vertically  integrated  advecticn-diffusion 
equation  for  a  salinity  or  water  quality  parameter  has  been  the  object 
of  a  few  investigations.  For  example,  Leendertse  and  Gritton  (1971) 
hindcasted  water  quality  parameters  in  Jamaica  Bay,  N.  Y. ,  and  Hess 
and  White  (1977)  modeled  the  dispersal  of  a  marked  fluid  in  Narragansett 
Bay  from  point  sources  within  the  bay. 

Three-Dimensional  Models 


16.  Three-dimensional  models  are  very  complex  and  only  a  few  have 
been  developed  to  a  point  where  simulations  are  feasible.  The  primary 
ones  are  the  model  investigations  by  Leendertse  and  Liu  ( 197 5 »  1977). 

This  model  uses  explicit  finite  difference  methods  and  consequently  is 
very  costly  in  computer  time  for  any  lengthly  simulation.  Horizontal 
grid  dimensions  are  similar;  therefore,  narrow  estuaries  are  difficult  to 
model.  Wang  (198c)  briefly  mentions  the  use  of  simi-implicit  methods  and 
mode  splitting  as  aids  to  efficiency  for  three-dimensional  models. 
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Methods  using  vertical  modes  instead  of  a  vertical  grid  have  been  devel¬ 
oped  by  Heaps  (197*0  for  shallow  coastal  seas.  This  model  is  linear  and 
does  not  include  the  advection  diffusion  equation  for  salinity  or 
density. 

17.  Caponi  (1976)  developed  a  complex  time-dependent,  three- 
dimensional  model  which  used  the  complete  three-dimensional  Navier- 
Stokes  equations  instead  of  the  Boussinesq  approximation. .  It  was  applied 
to  Chesapeake  Bay,  but  the  relatively  coarse  grid  used  seems  to  be  re¬ 
sponsible  for  some  anomalous  circulations. 

Box  Models 

18.  Officer  (1980)  has  critically  reviewed  box  models.  The 
estuary  is  divided  up  into  boxes  or  volumes  of  water  and  the  fluxes 
between  the  boxes  calculated  from  the  equations  of  conservation  of 
volume  and  salt.  The  resulting  calculation  is  very  simple  and  results 
in  a  steady-state  distribution  of  fluxes.  Elliott  (1975)  extended  the 
box  model  by  using  an  approximation  to  the  horizontal  momentum  equation 
to  calculate  velocities.  These  kinds  of  models  are  sometimes  used  as  a 
basis  for  water  quality  models  (Barrett  and  Mollowney  1972).  Numerical 
procedures  are  usually  based  on  Pritchard  (1969).  Officer  gives  examples 
where  two-dimensional  box  models  give  reasonable  results  when  compared 
with  two  studies  using  steady-state  numerical  models.  The  investiga¬ 
tions  are  the  uptake  of  silica  in  San  Francisco  Bay  (Peterson,  Festa, 
and  Conomos  1978)  and  the  theoretical  study  of  the  turbidity  maximum  in 
an  idealized  estuary  by  Festa  and  Hansen  (1978).  A  box  model  has  also 
been  used  to  assist  the  interpretation  of  nutrient  distributions  in 
Chesapeake  Bay  (Taft,  Elliott,  and  Taylor  1978). 

Summary 

19.  Table  1  summarizes  important  examples  of  the  different  kinds 
of  models  discussed  above.  It  is  not  an  exhaustive  list  of  all  the 
models  known  to  the  authors,  but  rather  shows  the  variety  of  approaches 
to  the  numerical  modeling  of  estuaries. 
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PART  III:  SEDIMENT  TRANSPORT  MODELS 


20.  Sediment  transport  and  the  interaction  of  sediment  particles 
and  turbulent  estuarine  flow  is  one  of  the  most  difficult  problems  in 
fluid  mechanics.  It  is  not  surprising  that  there  are  relatively  few 
sediment  transport  or  turbidity  models  compared  with  hydrodynamic  cir¬ 
culation  models.  Of  course,  a  sediment  transport  model  requires  a 
hydrodynamic  model  as  a  basis.  The  greater  complexity  in  the  hydro¬ 
dynamics  over  water  quality  models  arises  from  the  bottom  boundary  layer 
dynamics  where  important  processes  of  deposition  and  resuspension  take 
place . 

21.  Suspended  sediment  is  important  for  water  quality  models  be¬ 
cause  it  limits  light  penetration  and  thus  photosynthesis.  Also,  many 
pollutants  attach  readily  to  sediment  particles  and,  if  deposition  is 
taking  place,  may  be  removed  from  the  water  column.  However,  a  large 
storm  may  resuspend  large  quantities  of  polluted  sediments  which  could 
severely  affect  the  water  quality  of  the  estuary. 

22.  Festa  and  Hansen  (1978)  used  a  steady-state,  two-dimensional 
(longitudinal  and  depth)  model  for  the  circulation  in  an  idealized  es¬ 
tuary.  They  investigated  the  effect  of  settling  velocities  on  the  posi¬ 
tion  of  the  turbidity  maximum  assuming  no  flux  of  sediment  into  the  bed. 

23.  Odd  and  Owen  (1972)  have  developed  a  two-layer  model  of  the 
bed  flow  of  fluid  mud  in  the  Thames  Estuary  (United  Kingdom).  This 
model  is  time  dependent  and  solves  equations  for  bulk  flow  in  each  layer. 

2L.  Ariathurai  and  Krone  (1976a)  developed  a  two-dimensional 
(horizontal)  finite  element  model  for  the  depth- averaged  advection  diffu¬ 
sion  equation.  This  model  allows  for  deposition  and  resuspension;  the 
bed  is  modeled  as  a  number  of  thin  layers  ,  so  consolidation  of  the  bed 
can  be  accounted  for  (Ariathurai  and  Krone  1976b).  The  model  has  since 
been  developed  and  used  in  breadth-averaged  form  so  that  the  two  dimen¬ 
sions  are  longitudinal  and  vertical.  This  model  has  been  applied  to  the 
Savannah  Estuary  (Ariathurai,  MacArthur,  and  Krone  1977)* 

25.  Kuo,  Nichols,  and  Lewis  (1978)  used  a  two-dimensional  breadth- 
averaged  hydrodynamic  model  similar  to  that  of  Blumberg  (1975).  The 
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deposition  and  resuspension  are  modeled  by  simple  source  and  sink  terms 
in  the  bottom  layer  of  the  model.  The  model  is  qualitatively  compared 
with  data  from  the  Rappahannock  Estuary. 
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PART  IV:  WATER  QUALITY  MODELS 


Introduction 


2 6.  In  Part  II,  hydrodynamic  models  of  estuaries  were  discussed. 
These  models  are  the  basis  of  all  water  quality  models  which  attempt  to 
describe  the  distribution  of  constituents  in  an  estuary.  The  majority 
of  these  models  are  based  on  the  mass  balance  equations  of  continuity, 
momentum,  and  constitutent  conversation,  the  later  being  affected  by 
sources  p.nd  sinks  due  to  chemical  and  biological  interactions.  It  is 
these  source  and  sink  terms  that  are  the  main  subject  of  this  section. 

27.  Nearly  all  water  quality  models  are  one  dimensional  and  so 
use  the  const itutient  equation: 


3(AC.) 

_ 


3  (AUC .  ) 

- L_  + 

3x  3x 


(1) 


where 


A  =  cross-sectional  area 
Ch  =  concentration  of  the  ith  constituent 
t,x  =  time  and  the  longitudinal  dimension,  respectively 
U  =  cross-sectional  mean  velocity 

=  longitudinal  dispersion  coefficient 
S. .  =  source  and  sinks  of  constituent  i 

ij 

The  terms  U  and  A  are  calculated  from  a  one-dimensional  tidal  or 
steady-state  model,  and  E^  is  empirically  determined  by  matching  ob¬ 
served  longitudinal  distributions  of  salinity  (S.,  =  0)  using  Equation  1. 

^■0 


Solution  Techniques 


28.  Methods  for  the  solution  of  Equation  1  usually  involve  finite 
differences.  If  the  advection-diffusion  equation  is  solved  for  a  con¬ 
servative  substance  (S^  =  0),  then  stability  criteria  on  the  time  step 
depend  on  U  and  E^  .  The  maximum  allowable  time  step  is  often  much 
greater  than  the  usual  explicit  methods  for  solving  the  hydrodynamic 


equations.  Therefore,  if  Equation  1  is  solved  along  with  the  hydro- 
dynamic  equations  for  U  and  A  ,  it  may  not  be  necessary  to  calculate 
Equation  1  every  time  step.  Alternatively,  the  model  employing  a  set 
of  equations  may  be  solved  independently  of  the  hydrodynamics  where  U 
and  A  are  the  appropriately  time-averaged  output  from  the  hydrodynamic 
model.  In  some  cases,  the  hydrodynamic  model  is  steady-state  and  the 
water  quality  model  time  dependent.  Evidently  U  and  A  are  then  just 
constant  coefficients  in  Equation  1. 

29.  When  the  two-dimensional  (or  possibly  three-dimensional) 
version  of  Equation  1  is  being  employed,  then  spatial  averaging  of  the 
hydrodynamic  model  output  may  be  necessary  if  the  water  quality  model 
uses  large-scale  volume  elements.  For  example,  in  Kremer  and  Nixon’s 
(1978)  ecosystem  model,  initial  studies  of  the  hydrodynamics  in 
Narragansett  Bay  were  performed  using  the  two-dimensional  tidel  model 
of  Hess  and  White  (197M.  This  hydrodynamic  model  uses  32U  elements  of 
0.5  square  nautical  mile  (1.72  sq  km).  The  ecosystem  model  uses  eight 
different-sized  elements  to  represent  the  bay.  The  hydrodynamic  model 
was  used  to  study  the  exchange  of  material  between  the  eight  elements 
as  a  function  of  tide  height,  and  linear  regression  equations  were 
derived,  which  expressed  the  mass  exchange  between  elements  on  a  daily 
basis  as  a  function  of  tide  height.  Thus,  a  years  simulation  of  the 
exchanges  could  be  run  extremely  efficiently  compared  with  the  full 
hydrodynamic  model.  The  hydrodynamics  were  further  refined  by  using  a 
simple  box  model  calculation  (Pritchard  1969)  to  represent  two- layer 
exchanges  due  to  the  gravitational  circulation  in  the  bay.  Therefore, 
the  hydrodynamics  in  the  ecosystem  model  are  made  simple  and  efficient 
but  are  based  solidly  on  investigations  with  a  complex  time-dependent, 
two-dimensional,  vertically  integrated  hydrodynamic  model. 

30.  When  a  complex  set  of  source  and  sink  terms  are  part  of  a 
set  of  equations  describing  a  number  of  substances ,  then  the  stability 
of  the  set  of  the  differential  equations  has  to  be  analyzed  before  a 
solution  technique  for  the  source  and  sink  terms  can  be  determined. 

Since  usually  involves  no  spatial  or  time  gradients  said  is  usually 

linear  or  quasilinear,  the  spatial  dependence  of  the  left  hand  side  (LHS) 
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of  Equation  1  can  be  neglected  and  the  stability  or  stiffness  of  a  set 
of  ordinary  first-order  differential  equations  can  be  determined.  A  set 
of  stiff  ordinary  differential  equations  may  require  complex  iterative 
techniques  for  successful  solution.  The  more  simple  the  technique  for 
the  source  ard  sink  terms,  the  more  compa+ible  it  will  be  with  finite 
difference  techniques  used  for  the  LHS  of  Equation  1. 

31.  There  is  virtually  no  essential  difference  between  finite 
element  and  finite  difference  techniques  for  one-dimensional  problems. 
When  two  or  more  space  dimensions  are  involved,  the  two  types  of  tech¬ 
niques  differ  in  their  formulations.  Finite  difference  techniques  are 
more  common,  but  finite  element  techniques  have  been  used  for  the  two- 
dimensional  advection-dif fusion  equation  (Ariathurai  and  Krone  1976a). 

32.  In  some  water  quality  studies  a  more  sophisticated  hydro¬ 


dynamic  model  has  been  used  as  the  base  for  the  water  quality  calcula¬ 
tions  (Peterson,  Festa,  and  Conomos  1978;  Zitta,  Shindala,  and  Corey 
1977;  Winter,  Banse,  and  Anderson  1975;  Leendertse  and  Gritton  1971; 
Pitbaldo  and  Prince  1977).  For  example,  Zitta,  Shindala,  and  Corey 
(1977)  discuss  a  two-dimensional  vertically  well-mixed  steady-state 
estuary  model  for  describing  the  distribution  of  salinity,  nitrogen, 
biochemical  oxygen  demand  (BOD),  and  dissolved  oxygen  (DO)  concentrations 
in  a  system  of  interconnected  tributaries  and  bayous.  The  model  was 
used  to  predict  the  effect  of  alternative  wastewater  management  methods 
on  the  water  quality  of  the  estuary. 

Water  Quality  and  Ecosystem  Models 

33.  Historically,  the  analysis  of  water  quality  in  streams,  estu¬ 
aries,  and  lakes  has  concentrated  on  the  DO  and  the  BOD  due  to  waste 
loads.  The  BOD  is  primarily  responsible  for  the  consumption  of  oxygen 
in  an  estuary.  However,  there  are  many  other  interactions  that  affect 
DO,  which  tend  to  make  its  distribution  in  the  water  rather  complex. 
Against  this,  DO  is  a  good  index  of  the  general  health  of  the  aquatic 
environment.  A  review  of  oxygen  requirements  of  fish  is  given  by 
Alabaster  (1973)  and  Vernberg  (1971). 


34.  Organic  processes  in  estuaries  are  very  complex,  including 


primary  production,  respiration,  and  mineralization.  These  processes 
can  cause  significant  changes  to  DO  and  carbon  dioxide  concentration 
and  thus  alter  the  speciation  of  many  chemical  constituents.  An  estuary 
is  a  complex  interface  of  three  quite  separate  and  contrasting  environ¬ 
ments:  seawater,  fresh  water,  and  land.  The  estuary  is  very  dynamic 
and  the  turnover  rates  of  material  can  be  high.  These  processes  are 
dependent  on  fluctuations  of  the  marine,  limnetic,  and  terrestrial  in¬ 
fluences.  A  good  review  of  organic  processes  in  estuaries  is  given 
by  Head  (1976). 

35-  Perhaps  the  most  complete  view  of  an  aquatic  system  centered 
on  DO-BOD  is  given  by  Beck  (1978)  which  is  reproduced  in  Figure  1.  Of 
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Figure  1.  Factors  affecting  DO-BOD  interaction  (from 
Beck  1978,  courtesy  of  John  Wiley  &  Sons,  Ltd.  "Model¬ 
ing  of  Dissolved  Oxygen  in  a  Non-Tidal  Stream,"  by  M.  B. 
Beck  (A.  James,  ed.),  1978,  Wiley-Interscience,  New  York 


course,  the  degree  of  emphasis  on  the  various  constituents  and  the  inter¬ 
actions  included  in  the  model  vary  among  authors.  The  block  diagram  used 
to  specify  a  zooplankton,  phytoplankton,  nutrients,  and  DO-BOD  system 
is  given  by  DiToro  et  al.  (1977)  and  is  reproduced  in  Figure  2.  The 


Figure  2.  Kinetic  interactions  of  the  variables  of  an  estuarine 
water  quality  model  (from  DiToro  et  al.  1977*  courtesy  of  John 

Wiley  &  Sons,  Ltd.) 

nitrogen  cycle  is  specified  in  more  detail  and  the  DO-BOD-nitrogen 
interactions  is  in  less  detail.  The  water  quality  models  developed  by 
the  Water  Quality  Institute  in  Denmark,  which  are  reviewed  by  Dahl- 
Madsen  (1978),  are  also  complex  in  that  zooplankton,  phytoplankton, 
detritus,  nutrients,  and  DO  are  state  variables.  Processes  such  as 
mineralization  of  detritus  and  oxygen  consumption  by  the  sediment  are 
included.  The  models  also  keep  track  of  the  uptake  of  nitrogen  and 
phosphorus  in  the  phytoplankton  as  well  as  total  organic  carbon.  Fven 
though  the  DO-BOD  cycle  is  included,  the  main  interest  in  both  DiToro 
et  al.  (1977)  and  Dahl-Madsen  (1978)  is  in  the  problem  of  eutrophication 
due  to  excessive  algal  blooms  in  nutrient-rich  waters.  Another  complex 
ecosystem  model  which  includes  DO-BOD  is  by  Chen  and  Orlob  (1972,  1975); 
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this  model  also  includes  some  of  the  higher  trophic  levels  such  as  fish 
and  benthic  animals.  The  model  was  applied  to  eutrophication  problems 
in  Lake  Washington  and  San  Francisco  Bay.  These  models,  DiToro  et  al. 
(1977),  Dahl-Madsen  (1978),  and  Chen  and  Orlob  (1972),  are  effective 
ecosystem  models  which  have  been  developed  for  water  quality  problems 
concerning  DO  and  eutrophication.  However,  ecological  models  have  also 
been  developed  as  tools  for  investigating  processes  in  the  aquatic  eco¬ 
system  per  se.  For  example,  Kremer  and  Nixon  (1978),  in  their  ecosystem 
model  of  Narragansett  Bay,  have  a  very  complex  zooplankton  compartment 
and  consider  predation  by  the  benthos;  however,  they  do  not  use  DO-BOD 
as  part  of  the  system.  These  differences  reflect  the  goals  and  emphasis 
of  the  various  models  and  the  fact  that  there  is  no  universal  water 
quality  or  ecosystem  model  of  an  estuary. 

36.  The  majority  of  recent  sophisticated  water  quality  models  are 
based  on  three  fundamental  principles.  The  first  and  most  fundamental 
is  the  conservation  of  mass  as  given  by  Equation  1.  The  second  is  the 
kinetic  principle  in  which  generally  the  rate  of  change  of  constituent 
is  equal  to  the  product  of  constituents  that  affect  that  reaction  and  a 
rate  coefficient.  Thus,  the  rate  of  utilization  of  oxygen  by  bacteria 
is  linearly  proportional  to  the  BOD  load.  The  third  principle  is  the 
stoichiometry  of  chemical  transformations.  Whenever  the  conservation  of 
more  than  one  chemical  species  is  considered  in  a  model,  then  the 
stoichiometric  relationships  between  the  species  during  chemical  trans¬ 
formations  are  specified.  For  example,  if  along  with  the  conservation 
of  oxygen,  the  conservation  of  carbon  dioxide  was  modeled,  then  the 
ratio  of  oxygen  consumed  to  carbon  dioxide  generated  by  the  bacteria 
must  be  specified.  In  many  models  constant  stoichiometric  relationships 
are  implicitly  assumed  (i.e.,  DiToro  et  al.  1977).  However,  when 
phenomena  such  as  competitive  uptake  of  nutrients  by  algae  are  to  be 
modeled,  then  a  variable  stoichiometry  is  necessary. 

Water  quality  models  based  on  DO-BOD 

37-  The  majority  of  BOD- DO  models  have  evolved  from  the  study 
of  Streeter  and  Phelps  (1925 )•  They  assumed  that  the  balance  between 
DO  and  BOD  concentrations  is  the  result  of  two  processes  only:  the 


reaeration  of  the  stream  and  the  consumption  of  DO  in  the  oxidation  or 
decay  of  BOD  by  first-order  linear  kinetics.  This  model  still  has 
utility,  particularly  in  nontidal  streams  (Grantham,  Schaake,  and  Pyatt 
1971) •  Improvements  in  this  scheme  were  made  by  Dobbins  (196U)  and 
Dresnack  and  Dobbins  (1968)  who  added  source  and  sink  terms  for  DO  due 
to  plant  photosynthesis  and  respiration  or  decomposition  of  bottom  muds, 
and  sink  terms  for  BOD  due  to  sedimentation  and/or  adsorbtion.  These 
early  models  did  not  use  the  full  time- dependent  mass  balance  equation 
(Equation  l)  and  subsequent  work  corrected  these  deficiences  (Whitehead 
and  Young  1975;  Beck  1973;  Grantham,  Schaake,  and  Pyatt  1971;  O'Connor 
1967 ) •  A  further  step  in  the  development  of  DO-BOD  models  was  made  by 
O'Connor  (1967)  and  O'Connor  and  DiToro  (1970)  where  BOD  and  DO  utiliza¬ 
tion  were  divided  into  nitrogenous  and  carbonaceous  forms  and  rates.  The 
two  forms  of  BOD  are  associated  with  different  types  of  sewage  loading 
and  different  types  of  bacteria.  The  nitrogenous  form  is  often  asso¬ 
ciated  with  bottom  sediments  and  benthic  bacteria. 

38.  Further  progress  in  model  structure  has  been  made  by  the 
introduction  of  nonlinear  Monod  (19*+2)  or  Michaelis-Menten  kinetics. 

This  hyperbolic  growth  rate  has  been  used  as  a  description  of  algae 
growth  (Chen  1970;  Steele  197*+;  Thomann,  DiToro,  and  O'Connor  197*+; 

Kremer  and  Nixon  1978)  and  for  the  growth  of  bacteria  and  protozoa  in 
river  sediments  (Rutherford  and  O'Sullivan,  197*+)-  One  of  the  most 
comprehensive  estuarine  water  quality  models  is  that  by  Thomann,  DiToro, 
and  O'Connor  (197*+)  and  later  developments  by  DiToro  et  al.  (1977) • 

In  these  models  algae  is  a  major  compartment  affecting  both  the  nitrogen 
cycle  and  BOD-DO.  It  is  noted  that  in  all  these  models  the  rate  con¬ 
stants  of  decay  and  growth,  whether  from  linear  kinetics  or  nonlinear 
Monod  kinetics,  are  highly  sensitive  to  temperature.  Additionally, 
these  models  require  forcings  to  be  specified  as  functions  of  time. 
Forcings  include  temperature,  light  (for  photosynthesis),  wind  (for 
reaeration),  waste  loading  (for  BOD),  as  well  as  the  time-dependent 
tidal  and  nontidal  estuarine  flow. 

39.  Monod  hyperbolic  functions  which  relate  growth  to  algae  or 
bacteria  to  the  substrate  are  almost  universally  used  in  water  quality 
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and  ecologic  models.  However,  among  biologists  there  are  some  doubts 
whether  Monod  kinetics  are  correct,  particularly  in  multinutrient  en¬ 
vironments  with  hetrogeneous  algal  populations  (Platt,  Denman,  and 
Jassby  1977)*  Models  have  been  developed  by  Bierman  et  al.  (197M  for 
phytoplankton  species  competition  in  a  lake.  Cell  growth  was  treated 
as  a  two-step  mechanism  that  involves  separate  nutrient  uptake  and 
cell  synthesis  processes  and  includes  an  intermediate  storage  capability. 
The  mode  of  cell  growth,  i.e.,  assimilation  and  subsequent  ingestion  of 
nutrients,  was  thought  to  be  more  important  in  modeling  growth  dynamics 
than  are  the  extremely  complex  chemical  conversions  within  the  cell. 

UO.  It  appears  the  more  complex  the  model  of  DO-BOD,  the  more 
kinetic  rate  coefficients  need  to  be  determined;  however,  it  is  more 
likely  that  these  rates  can  be  determined  independently  through  labora¬ 
tory  measurements  of  that  process.  For  example,  if  algal  respiration  is 
included  in  the  model  structure,  then  measurements  in  the  laboratory  or 
in  situ  may  be  able  to  determine  the  range  of  values  for  the  rate  co¬ 
efficient  for  which  the  model  must  produce  reasonable  results.  If  a 
number  of  processes  are  lumped  into  one  coefficient  as  in  the  Streeter 
and  Phelps  (1925)  model,  then  it  is  difficult  to  determine  the  decay- 
coefficients  independently  of  the  simulation. 

Ul.  This  point  leads  to  the  whole  question  of  simulation  and 
verification  of  water  quality  models.  In  almost  all  cases,  it  is  nec¬ 
essary  to  determine  some  coefficients  by  fitting  the  model  output  to 
observations  in  the  estuary  or  river.  The  determination  of  the  longitud¬ 
inal  dispersion  coefficient  Ex  has  already  been  mentioned.  This  pro¬ 
cess  is  known  as  calibration.  The  hydrodynamic  part  and  the  kinetic 
parts  of  the  model  are  regarded  as  being  independent  of  each  other. 

Once  a  model  has  been  calibrated,  then  it  must  be  verified  against  a 
completely  independent  set  of  data  with  no  change  in  the  parameters. 

Only  then  can  the  model  be  reliably  used  in  a  forecast  mode.  For  exam¬ 
ple,  DiToro  et  al.  (1977)  calibrated  their  phytoplankton  model  for  the 
Potomac  Estuary  by  using  1968  data.  They  verified  it  against  1969  data 
making  only  the  required  changes  in  the  forcing  functions,  riverflows, 
light,  waste-loading,  etc.  They  achieved  satisfactory  agreement  which 
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allowed  confidence  that  the  major  physical,  chemical,  and  biological 
processes  were  satisfactorily  included  in  the  model  structure. 

b2.  The  models  discussed  so  far  are  deterministic  and  parameterize 
the  internal  structure  of  the  system.  The  alternative  approach  is  the 
black  box  approach  which  is  often  used  in  rivers  where  time  series  of 
water  quality  variables  measured  at  an  upstream  site  are  used  to  pre¬ 
dict  variables  at  a  downstream  site.  This  is  often  useful  when  combined 
with  system  analysis  as  advocated  by  Thomann  (1972).  It  allows  the  use 
of  various  modeling  hypotheses,  within  the  black  box,  to  be  applied  to 
a  reach  of  a  river  and  evaluated  against  data.  A  case  study  on  the 
river  Cam  in  England  is  given  by  Beck  (1978).  This  stochastic  modeling 
approach  would  be  difficult  to  apply  to  an  estuary  with  its  reversing 
tidal  flows. 

^3.  A  further  question  that  arises  in  the  use  of  models  for 
managing  aquatic  systems  concerns  time  dependency.  Often  steady-state 
models  are  used  as  a  forecast  of  the  long-term  state  the  estuary  or 
river  will  reach  under  various  loadings.  The  time-dependent  models  are 
often  used  to  study  more  short-term  variability  of  the  estuary.  Of 
course,  the  latter  require  a  great  deal  more  computer  resources  than  the 
former,  and,  therefore,  the  former  are  often  a  cheap  way  of  investigating 
many  different  scenarios.  There  is  a  trend  at  present  towards  using 
time-dependent  models  for  all  situations,  particularly  for  estuaries 
which  are  highly  dynamic.  Harlemann  (1977)  has  pointed  out  that  use  of 
steady-state  models  of  complex  systems  can  lead  to  serious  errors  when 
compared  to  time-dependent  models. 

kk.  Reviews  of  water  quality  modeling  may  be  found  in  the  books 
edited  by  James  (1978);  Canale  (1976);  Biswas  (1976);  and  Hatzinger, 
Vanlelgveld,  and  Zoeteman  (1978);  and  in  review  articles  by  DiToro  et 
al.  (1977);  Hahn  and  Schreiner  (1978);  and  O'Connor,  Thomann,  and  DiToro 
(1977).  Perhaps  one  of  the  most  successful  uses  of  a  simple  kinetic 
model  to  assist  in  the  clean  up  of  an  estuary  is  the  study  of  the  Water 
Pollution  Research  Laboratory  in  England  on  the  Thames.  This  is  assessed 
by  Barrett,  Mollowney,  and  Casapieri  (1978). 


Water  quality  models 
not  involving  DO-BOD 

45.  Some  modeling  studies  have  been  performed  on  estuaries 
where  the  primary  interest  is  not  DO-BOD  or  the  higher  trophic  levels, 
but  the  distribution  of  silica.  Silica  is  utilized  by  algae  in  the 
estuary  and  constant  uptake  rates  are  assumed  which  may  vary  with  depth. 
A  more  complex  depth-dependent,  steady-state,  two-dimensional  circula¬ 
tion  model  is  used  (Festa  and  Hansen  1976),  and  the  relations  between 
advection,  mixing,  and  uptake  are  more  completely  characterized  (Peter¬ 
son,  Festa,  and  Conomos  1978).  Recently,  Rattray  and  Officer  (1979) 
showed  that  very  similar  results  could  be  obtained  for  San  Francisco 
Bay  by  using  a  variation  on  the  box  model  approach  to  characterize  the 
circulation  and  mixing.  This  approach  is  considerably  less  costly  than 
using  the  full  numerical  model.  Silica  uptake  has  also  been  studied  for 
the  Scheldt  Estuary  (The  Netherlands)  by  Wollast  (1977)  using  a  simple 
model . 

46.  An  interesting  approach  to  water  quality  modeling  in  a  highly 
polluted  estuary,  where  nutrient  limitation  is  not  a  factor,  is  taken 

by  Nihoul  et  al.  (1979)  and  Billen  and  Smitz  (1977).  Each  chemical  is 
considered  in  oxidized  and  reduced  form  through  reversible  reactions. 
Thus,  the  two  forms  together  are  conservative.  The  estuary  is  con¬ 
sidered  to  be  in  thermodynamic  equilibrium  and  redox  potentials  are 
involved  as  part  of  the  solution.  The  final  relationship  between  the 
eight  separate  reactions  is  made  by  considering  the  global  oxygen  budget. 
This  makes  the  solution  technique  particularly  elegant  since  the  redox 
potentials  which  control  the  reactions  are  known  from  laboratory  experi¬ 
ments.  The  authors  admit,  however,  that  thermodynamic  equilibrium  is 
not  necessarily  fulfilled  throughout  the  estuary. 

47.  In  contrast  to  models  of  BOD  and  nitrogen  loading  of  estu¬ 
aries,  there  have  been  few  models  of  the  fate  and  effects  of  heavy  metals 
One  such  model  is  that  by  Jorgensen  (1979)  which  is  based  on  the  concept 
of  trophic  length  levels  recently  introduced  by  Thomann  (1979)  where  a 
length  scale  is  introduced  into  the  equations  as  an  independent  variable 
which  represents  the  length  of  time  a  substance  spends  in  a  particular 
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trophic  level.  The  fate  of  metals  in  an  estuarine  environment  is  very 
complicated,  involving  precipitation  of  metals  by  seawater,  adsorbtion 
by  suspended  sediment,  and  biological  uptake  by  algae  and  benthic 
organisms.  The  geochemistry  of  estuarine  trace  metals  is  discussed  by 
Turekian  (1977)-  A  study,  using  a  model  based  on  population  dynamics, 
of  the  uptake  of  heavy  metals  by  benthic  algae  has  been  performed  by 
Seip  (1979)  and  Seip  et  al.  (1979)*  This  study  compares  the  differences 
found  between  a  polluted  and  nonpolluted  Norwegian  fjords. 

Time  and  Length  Scales 


U8.  All  water  quality  and  ecologic  models  are  forced  systems. 

An  important  part  of  the  construction  of  models  is  the  question  of 
resolution  of  the  important  time  and  length  scales.  To  a  large  degree, 
the  system  will  be  governed  by  the  scales  of  the  hydrodynamics.  Impor¬ 
tant  scales  range  from  the  turbulent,  on  the  order  of  a  few  centimetres 
and  a  few  seconds;  tidal,  12  to  25  hr  and  a  few  kilometres;  wind,  a  few 
days  and  a  few  tens  of  kilometres;  to  riverflow,  one  to  several  months 
(flushing  times)  and  tens  to  hundreds  of  kilometres.  Of  course,  these 
scales  are  functions  of  the  characteristic  geometry  (longitudinal,  width, 
and  depth)  of  the  estuary.  Modern  time-dependent  hydrodynamic  models 
normally  have  time  steps  on  the  order  of  a  few  minutes  to  1  hr  and 
so  should  be  able  to  resolve  all  the  above  scales  except  the  turbulent, 
which  is  parameterized  by  eddy  viscosity  and  diffusion  coefficients. 

The  degree  to  which  they  do  so  depends  on  the  sophistication  of  the  model 
and  the  number  of  space  dimensions  considered. 

U9.  Biological  and  chemical  processes  also  have  characteristic 
length  and  time  scales,  some  of  which  may  be  different  from  the  physical 
scales.  The  main  additional  forcing  mechanism  is  that  of  light  which 
has  diurnal  and  seasonal  variability  and  also  a  depth  scale  associated 
with  extinction  by  the  water.  If  point  and  nonpoint  sources  for  mate¬ 
rial  (such  as  BOD  waste)  are  separate  from  the  freshwater  flow,  then 
their  characteristic  time  and  length  scales  also  need  to  be  considered. 
Ford  and  Thornton  (1979)  have  reviewed  time  and  length  scales  as  they 
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affect  ecological  modeling  of  lakes.  Patchiness  of  plankton  in  the 
marine  environment  has  received  some  study  because  the  nature  of  the 
variance  is  important  for  predators  such  as  fish  larvae.  An  example  is 
the  investigation  by  Steele  and  Henderson  (1977)  which  confirms  that 
important  scales  appear  to  be  on  the  order  of  kilometres  and  days  for 
the  lower  trophic  levels  in  the  North  Sea.  At  horizontal  scales  less 
than  100  m,  chlorophyll  spectra  were  similar  to  the  turbulent  spectra 
of  temperature.  The  1-day  time  step  is  very  common  for  ecological 
models  (DiToro  et  al.  1977;  Chen  and  Orlob  1972);  however,  Jamart  et 
al.  (1977)  show  that,  for  a  depth-  and  time-dependent  model  of  phyto¬ 
plankton  growth  in  the  northeastern  Pacific,  due  to  the  nonlinear  nature 
of  the  light  forcing,  time  steps  of  6  min  are  required  during  times  of 
rapid  light  changes.  In  some  systems  (i.e.  Puget  Sound)  plankton  cells 
may  double  in  number  in  a  few  hours. 

50.  An  example  of  the  increasing  sophistication  of  ecological 
models  in  taking  account  of  the  time  scales  of  the  physical  medium  is 
given  by  Kemp  and  Mitsch  (1979).  They  suggest  tha,t  the  natural  time 
variability,  due  to  advection  and  turbulence,  of  the  euphotic  layer  of 

a  water  body  modifies  interactive  pressures  between  plankton  populations. 
The  growth  and  distribution  of  a  plankton  population  are  to  some  extent 
controlled  by  the  water  movements;  for  turbulent  frequencies,  which  are 
of  the  same  order  as  growth  rates,  heterogeneous  populations  are  shown 
to  coexist  in  the  model.  Without  the  turbulence,  one  population  would 
out  compete  the  others  and  completely  and  'unrealistically  dominate  the 
water  body. 

51.  The  chemical  and  biological  equations  of  water  quality  models 
have  their  own  characteristic  time  scales  determined  by  the  kinetic  rate 
coefficients.  These  time  scales  are  usually  on  the  order  of  1  to  10 
days;  however,  the  linear  kinetics  of  the  model  are  the  result  of  quasi- 
reversible  chemical  reactions  that  may  take  place  very  rapidly.  DiToro 
(1976)  has  given  some  theoretical  basis  by  which  relatively  slow  linear 
kinetics  are  satisfactory.  The  thermodynamic  equilibrium  argument  used 
by  Nihoul  et  al.  (1979)  that  employs  redox  potentials  is  very  similar 

to  that  of  DiToro  (1976). 
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Summary 

52.  Table  2  summarizes  examples  of  different  kinds  of  models 
used  for  water  quality  in  estuaries  that  have  been  discussed  above. 
Again,  this  is  a  list  of  important  models  that  are  representative  of 
different  approaches  and  is  not  meant  to  be  exhaustive. 
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PART  V:  ESTUARINE  ECOSYSTEM  MODELS 


53.  In  the  previous  section,  a  number  of  water  quality  models, 
which  were  effectively  applied  estuarine  ecosystem  models,  were  dis¬ 
cussed.  In  this  section,  estuarine  ecosystem  models  are  surveyed  which 
attempt  to  further  the  knowledge  of  the  interactions  between  the  lower 
trophic  levels  and  the  various  forcing  mechanisms.  A  major  difference 
between  these  models  and  water  quality  models  is  their  historical  line 
of  development  from  marine  ecologic  studies. 

Marine  Models 


5*4.  Modeling  primary  production  in  the  sea  has  a  long  history, 
which  has  been  reviewed  by  Cushing  (1975)  in  his  book  on  marine  ecology 
and  fisheries.  He  stresses  the  abandonment  of  the  agricultural  model 
based  on  limiting  nutrients  for  the  predator-prey  model  and  the  regenera¬ 
tion  of  nutrients  through  rapid  recycling  of  detritus.  The  majority  of 
marine  ecosystem  models  have  concentrated  on  the  nutrient-phytoplanktcn- 
zooplanktcn  interactions  where  often  there  is  only  one  compartment  for 
each  trophic  level,  but  transfers  between  compartments  are  complex 
(Steele  197*0.  The  reason  for  this  is  that  the  upper  layers  of  the  ocean 
are  relatively  homogeneous  (compared  with  terrestrial  systems)  and  phyto¬ 
plankton  species  show  only  minor  differences  in  basic  chemical  composi¬ 
tion.  It  seems  reasonable  to  lump  all  phytoplankton  species  into  one 
compartment.  These  kinds  of  investigations  do  not  involve  any  explicit 
spatial  variation  of  state  variables;  however,  this  basic  approach, 
which  has  been  developed  from  the  pioneering  studies  of  Fleming  (1939) > 
Riley  (19 b6,  19*47 ) 5  and  Riley,  Stommel,  and  Bumpus  (19*49)  >  has  been 
adapted  to  allow  spatial  as  well  as  time  variability  in  systems  which 
have  large  physical  gradients  such  as  upwelling  regions  (Walsh  1971; 
O'Brien  and  Wroblewski  1973)  and  estuaries  (Winter,  Banse,  and  Anderson 
1975;  Kremer  and  Nixon  1978).  The  alternative  approach  advocated  by 
Patten  (1968)  uses  very  complex  model  structures  involving  all  critical 
components,  which  may  involve  many  species,  but  the  transfers  between 


components  are  generally  simple,  linear,  and  empirical.  This  kind  of 
model  tends  to  be  very  theoretical  because  data  are  not  available  for 
calibration  of  the  transfers  and  has,  therefore,  been  little  used  for 
aquatic  systems. 

55.  The  higher  trophic  levels,  such  as  fish  or  shellfish,  have 
been  very  rarely  coupled  to  ecosystem  models  of  the  lower  trophic  levels. 
Green  (1975)  was  among  the  first  to  model  an  entire  ecosystem  with  her 
model  of  the  Ross  Sea.  Models  of  population  dynamics  of  fisheries  are 
used  for  fishery  management  (Beverton  and  Holt  1957);  examples  of  this 
approach  applied  to  estuaries  are  given  by  Wyatt  and  Quasim  (1973)  and 
Seip  et  al.  (1979). 


Estuarine  Models 


56.  Models  of  estuarine  systems  have  tended  to  lag  behind  those 
of  the  marine  system.  Many  of  the  more  recent  modeling  efforts  in  the 
estuarine  system  have  been  theoretical  in  nature,  i.e.,  have  attempted 
to  further  the  understanding  of  nature.  One  of  the  earliest  estuarine 
ecological  models  was  that  of  Pomeroy  et  al.  (1972).  This  model  was 
designed  to  simulate  phosphorus  flux  in  several  Georgia  estuaries. 

Bahr  (197^)  developed  both  linear  and  nonlinear  models  of  energy  flux  in 
Georgia  oyster  beds.  A  similar  model  by  Dame,  Vernbert,  and  Bonnell 
(1977)  on  South  Carolina  oyster  beds  emphasized  sensitivity  analysis, 
i.e.,  critical  components  in  the  system.  Show  (1979)  developed  a 
hydrological-biological  model  to  explain  zooplankton  distributions  in  a 
small  lagoon  on  the  Texas  coast.  Ferguson  and  Adams  (1979)  simulated 
the  trophic  dynamics  of  epifauna  and  juvenile  fish  of  an  eelgrass  commu¬ 
nity  in  order  to  test  their  response  to  temperature. 

57.  Probably  the  two  most  important  ecosystem  models  of  the  lower 
trophic  levels  of  estuaries  are  by  Winter,  Banse,  and  Anderson  (1975) 
and  Kremer  and  Nixon  (1978).  The  ecosystems  are  contrasted  primarily  by 
their  physical  dimensions  and  hydrodynamics.  Winter,  Banse,  and  Anderson 
(1975)  constructed  a  model  of  primary  production  for  the  central  basin 

of  Puget  Sound,  which  is  a  deep  fjord  with  strong  surface  layer 


stratification.  The  Kremer  and  Nixon  (1978)  model  is  for  Narragansett 
Bay,  which  is  a  shallow,  wide,  and  partially  mixed  estuary.  It  is  note¬ 
worthy  that  both  models  are  based  on  extensive  investigations  of  the 
ecosystems  carried  out  by  many  different  investigators  over  many  years. 

58.  The  hydrodynamic  basis  of  Winter,  Banse,  and  Anderson  (1975) 
is  an  analytical,  steady-state,  longitudinal,  and  depth-dependent  model 
of  the  circulation  and  stratification  as  it  is  Influenced  by  tidal 
mixing,  river  discharge,  and  wind.  The  primary  production  model  was 
forced  by  observed  nitrate  and  zooplankton  concentration  and  solar  in¬ 
solation.  Phytoplankton  concentrations  were  simulated  for  April-June 
19 66  and  April-May  19&7-  It  was  found  that  estuarine  circulation  with 
its  associated  upwelling  from  depth  supplies  the  euphotic  zone  with  algal 
seed  stock  and  replenishes  exhausted  supplies  of  essential  nutrients  so 
that  nutrient  limitation  is  rarely  a  factor  in  algal  growth  rates.  It 
was  also  found  from  the  model  that  herbivore  grazing  and  algal  sinking 
were  of  secondary  importance.  The  investigation  is  a  good  example  of 
how  a  model  can  assist  the  quantitative  analysis  of  the  dynamics  of 
phytoplankton  blooms.  New  relationships  between  system  components,  such 
as  the  strong  influence  of  physical  environmental  factors  on  phytoplank¬ 
ton  growth,  are  often  perceived,  and  experimental  priorities  are 
classified. 

59-  The  state  of  the  art  in  estuarine  ecosystem  modeling  is 
probably  represented  by  the  Kremer  and  Nixon  (1978)  model  of  Narragansett 
Bay.  A  detailed  energy  flow  model  of  the  bay  was  developed  based  on 
information  taken  from  various  sources  over  a  20-year  period.  This  con¬ 
ceptual  model  was  much  more  simple  than  nature,  yet  too  complicated  to 
simulate.  The  model  eventually  chosen  for  simulation  emphasized  the 
plankton-nutrient  components  of  this  system  (Figures  3  and  *0.  The  six 
primary  state  variables  included  phytoplankton;  smaller,  primarily  herbi¬ 
vorous  zooplankton;  ammonia;  nitrate  and  nitrite;  phosphate;  and  sili¬ 
cate.  In  contrast  to  Winter,  Banse,  and  Anderson  (1975).  the  nutrient 
cycle  and  zooplankton  are  included  in  the  model  and  the  strong  inter¬ 
actions  and  feedback  between  the  three  major  compartments  are  very 
important  in  the  model  formulation.  Secondary  compartments  that  were 
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Figure  3.  A  complex,  but  still  greatly  simplified,  energy  flow  diagram  for  the 
Narragansett  Bay  ecosystem  on  a  summer  day.  This  conceptual  model  was  a  first 
step  in  the  process  of  abstraction  that  led  to  development  of  the  numerical 
model  (from  Kremer  and  Nixon  19T8»  courtesy  of  Springer-Verlag ,  Inc.) 
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Figure  4.  Final  energy  flow  diagram  and  the  conceptual 
framework  for  the  numerical  model  of  the  Narragansett 
Bay  ecosystem.  While  much  of  the  complexity  of  Fig¬ 
ure  3  has  been  omitted,  some  detail  has  been  added  in 
the  zooplankton  compartment.  Only  the  phytoplankton, 
zooplankton,  and  nutrient  compartments  are  fully  simu¬ 
lated  with  mechanistic  detail.  This  same  conceptual 
model  is  used  in  each  of  the  eight  elements  of  the  hay, 
with  appropriate  values  for  each  compartment  and  forc¬ 
ing  function  (from  Kremer  and  Nixon  1978,  courtesy  of 
Springer-Verlag,  Inc.) 

thought  to  he  important  influences  on  the  plankton-nutrient  system  were 
included  by  forcing  appropriate  patterns  when  necessary  to  express 
feedbacks  and  internal  controls.  The  primary  control  exerted  by  the 
higher  trophic  levels  was  predation  by  the  benthos.  Narragansett  Bay 
was  divided  into  eight  areas,  and  the  coupling  to  the  hydrodynamics  is 
described  in  Part  IV.  Simulations  done  in  each  area  account  for  rela¬ 
tively  large-scale  spatial  heterogeneity  within  the  bay. 

60.  The  model  had  various  limitations  that  were  discussed  by 
Kremer  and  Nixon  (1978).  The  model,  at  best,  represents  the  current 
understanding  of  the  Narragansett  Bay  system.  In  its  present  design. 
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the  model  cannot  change  to  different  system  states  that  may  result  from 
large  increases  in  inorganic  inputs,  drastic  changes  in  salinity,  etc. 
Hence,  the  model  is  limited  as  a  management  tool  to  those  responses  of 
the  system  due  to  small  changes  in  parameters  and  processes  that  are 
specifically  included  in  the  formulation  of  the  model.  However,  the 
model  is  essentially  a  complex  tool  designed  to  assist  in  the  analysis 
of  a  complex  estuarine  ecosystem  process  and  to  integrate  the  data  into 
an  overall  perspective  and  framework.  As  the  authors  point  out,  a  model, 
no  matter  how  realistic,  cannot  offer  escape  from  value  judgments  and 
long  links  of  supposition  for  the  decisionmakers.  An  important  point 
was  that,  if  a  management  tool  is  contemplated,  then  the  management 
objectives  and  related  questions  to  be  posed  to  the  model  should  be 
specified  from  the  outset. 

Estuarine  Management  Models 

6 1.  Empirical  models  are  often  used  in  management-related  studies 
and  were  discussed  in  Part  IV.  The  models  of  O'Connor,  DiToro  and 
Mancine  (1972);  O'Conner,  DiToro,  and  Thomann  (1975);  and  Chen  and  Orlob 
(1972,  1975)  have  been  applied  to  water  quality  problems  in  various 
aquatic  systems  around  the  country.  O'Connor  and  his  associates  over¬ 
come  the  lack  of  data  inherent  in  model  building  by  adapting  their  mod¬ 
els  to  the  existing  data  base,  whatever  it  may  be,  rather  than  the 
reverse  (Orlob  1975)-  This  tends  to  lead  to  aggregation  (i.e.,  temporal 
and  spatial  averaging);  however,  the  model  tends  to  be  amplified  con¬ 
sistent  with  the  data  base  and  the  level  of  detail  required  by  decision¬ 
makers.  Kelly  (1975)  and  Kelly  and  Spafford  (1977)  developed  their 
ecological  model  of  the  Delaware  River  estuary  specifically  to  answer 
management  questions  with  regard  to  water  quality  problems.  One  of  the 
goals  of  this  exercise  was  to  show  how  ecological  models  could  be 
incorporated  into  wider  scope  management  (optimization)  models.  The 
authors  concluded,  however,  that  a  lack  of  appropriate  data  is  the  major 
limiting  factor  in  constructing  and  validating  ecological  models  used 
for  management  purposes. 
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Fishery  Models 


62.  Another  type  of  practical  model  widely  used  in  the  past  is 
the  fish  population  model.  These  models  have  generally  emphasized  single 
species  and  are  not  true  ecological  models  in  the  sense  that  they  do  not 
show  relationships  between  species.  In  view  of  the  fact  that  most 
estuarine  models  have  neglected  the  higher  trophic  levels  except  in  a 
cursory  way,  applications  of  the  information  developed  from  these  mod¬ 
els  may  prove  useful  in  building  estuarine  ecosystems  models.  Schaaf 
(1975)  suggests  that  traditional  fishery  biology  will  benefit  from  a 
broader  ecosystem  viewpoint  and  that  extensive  data  and  population  model¬ 
ing  of  fishery  biology  can  play  an  important  role  in  refining  ecological 
models. 


PART  VI:  MARINE  WETLAND  ECOSYSTEM  MODELS 


63-  Marsh  ecosystems  are  biologically,  chemically,  and  physically 
complex  systems.  Figure  5  shows  one  view  of  the  relationships  that 
exist  between  these  parameters  in  terms  of  the  carbon  flux  through  the 
system.  Three  zones  of  ecological  interaction  are  depicted  in  this 
diagram:  (l)  terrestrial  interactions  in  the  abovewater  part  of  the 

marsh;  (2)  aquatic  interactions  within  the  water  in  the  tidal  creeks  and 
on  the  marsh  during  high  tide;  and  (3)  interactions  within  the  largely 
anaerobic  sediments  (Wiegert  and  Wetzel  1979)-  The  earliest  marsh 
models  tended  to  concentrate  on  interactions  within  just  one  of  these 
zones;  more  recent  modeling  efforts  have  incorporated  interactions  in 
all  three  of  the  zones,  although  at  various  levels  of  complexity. 


Figure  5.  Diagram  of  carbon  flow 
in  a  salt  marsh.  Carbon  is 
viewed  as  either  living  or  non¬ 
living  and  the  marsh  is  divided 
into  the  air,  water,  and  sediment 
zones  (courtesy  of  Weigert  and 
Wetzel  1979) 
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Marsh  Zonal  Models 


1 


I 

4 


6k.  Teal  (1962)  was  among  the  first  to  use  a  conceptual  model  to 
demonstrate  carbon  flow  through  a  system.  This  model  concentrated  on 
the  biological  components  in  the  emergent  portion  of  the  marsh.  The 
model  was  not  developed  beyond  the  conceptual  stage. 

65.  Studies  by  Nixon  and  Oviatt  (1973)  emphasized  energy  flow  in 
the  marsh  creeks  and  embayments  that  are  linked  with  the  tall  grass 
stands  and  that  serve  as  tidal  pathways  between  the  emergent  marsh  and 
the  estuaries  and  offshore  waters.  Conceptual  models  were  constructed 
for  these  habitats  for  both  the  stammer  and  winter.  A  simulation  model 
was  developed  which  provided  some  quantitative  detail  in  the  conceptual 
model  of  diurnal  dissolved  oxygen  in  the  tidal  creeks  and  embayments. 

66.  Although  the  oxygen  model  included  both  biological  and  physi- 

* 

cal  inputs,  the  amount  of  data  available  to  construct  the  model  and  the 
need  for  simplification  placed  limits  on  the  model  (Nixon  and  Oviatt 
1973).  The  model  was  used  to  simulate  the  effects  of  sewage  input  and 
temperature  increases  on  oxygen  concentrations  in  the  embayments.  Tidal 
inputs  were  found  to  have  marked  effects  on  the  oxygen  levels  in  the 
tidal  creeks  and  embayments.  The  model  in  itself  is  useful  in  various 
respects,  but,  because  of  the  limits  inherent  in  the  model,  could  best 
be  utilized  when  coupled  to  a  total  ecosystem  model.  Predictions  affect¬ 
ing  biotic  responses  in  a  system  could  then  be  based  on  model  results 
rather  than  extrapolating  to  the  biota  from  data  on  oxygen  concentrations 
in  the  system. 


Integrated  Marsh  Models 

67.  More  recent  models  have  become  more  complex  in  their  treat¬ 
ment  of  marsh  ecosystems.  Wiegert  and  Wetzel  (1979)  have  developed  a 
model  of  a  Georgia  Spartina  marsh  that  includes  seven  biotic  and  seven 
abiotic  compartments  (Figure  6).  This  model  represents  a  total  eco¬ 
system  approach  in  that  it  incorporates  biological,  chemical,  and 
physical  interactions  in  the  three  ecological  zones  discussed  earlier. 
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Figure  6.  Pathways  of  carbon  flux  between 
the  components  of  a  ll-compartment  salt 
marsh  model  continuing  the  structural  and 
functional  view  outlined  in  Figure  5. 

Solid  lines  designate  donor-controlled 
pathways.  Dotted  lines  designate  f-luxes 
that  respond  to  changes  in  either  or  both 
donor  and  recipient  (courtesy  of  Weigert 
and  Wetzel  1979) 

Besides  the  ll  compartments,  the  model  has  91  parameters  and  44  fluxes. 

68.  The  model  is  presently  in  its  third,  and  final,  version. 
Earlier  versions  of  the  model  identified  21  parameters  sensitive  enough 
to  change  the  final  system  steady-state  organic  accumulation  by  0.5  per¬ 
cent  (Wiegert  et  al.  1975).  These  sensitive  parameters  included  the 
components:  Spartina  shoots  and  roots,  algae,  particulate  organic  sedi¬ 

ments,  sediment  saprophages  (both  aerobic  and  anaerobic),  and  open-water 
heterotrophs.  These  earlier  models  also  suggested  that  the  scarcity  of 
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space  (as  opposed  to  scarcity  of  material  resources)  is  the  major  limit¬ 
ing  factor  for  most  components. 

69.  The  most  recent  version  of  the  model  produced  an  acceptable 
simulation  of  the  salt  marsh  system.  The  value  of  this  model  is  that 
it  has  attempted  to  incorporate  into  the  model  all  of  the  known  physical 
and  chemical  processes  which  influence  carbon  fluxes.  However,  in  order 
to  simplify  the  model,  compartments  had  to  be  generalized  (e.g.,  grazers). 
This  removed  some  of  the  usefulness  of  the  model  in  predicting  impacts 
to  specific  biotic  components  (e.g.  birds).  However,  the  model  has  been 
successful  in  satisfying  its  initial  objective  of  describing  and  explain¬ 
ing  the  pathways  and  dynamics  of  carbon  flow  in  the  marsh. 

TO.  Hopkinson  and  Day  (1977)  developed  individual  conceptual 
models  for  carbon  and  nitrogen  flows  in  a  Louisiana  salt  marsh.  These 
two  conceptual  models  were  then  combined  to  produce  a  single  simulation 
model  for  the  marsh  system.  This  simulation  model  incorporated  seven 
compartments  from  the  carbon  flow  model  and  three  nitrogen  compartments. 
These  latter  compartments  act  in  many  places  to  control  carbon  flows. 

As  in  the  Wiegert  and  Wetzel  model,  it  was  necessary  to  aggregate  com¬ 
ponents  of  the  ecosystem  into  single  compartments.  Hopkinson  and  Day 
(1977)  recognized  that  this  aggregation  affected  the  realism  of  the 
model,  which  was  evident  in  the  model  results.  Probably  the  primary 
value  of  this  modeling  exercise  is  that  it  suggested  new  areas  for  fur¬ 
ther  research  and  demonstrated  the  usefulness  of  models  as  a  research 
tool  besides  being  an  end  in  itself. 

71.  Reimold  (197*0  and  Zieman  and  Odum  (1978)  have  developed  mod¬ 
els  that  simulated  specific  components  of  the  system  rather  than  taking 
the  ecosystem  approaches  described  above.  Reimold 's  model  was  sensitive 
to  discrete  changes  in  Spartina  alterniflora.  The  model  consisted  of 
five  compartments  including  water,  sediments,  Spartina  detritus  and 
microflora,  and  detritus  feeders.  The  model's  usefulness  lies  in  the 
capability  of  the  model  to  predict  the  effects  of  perturbations  on  the 
Spartina  system. 

72.  Zieman  and  Odum's  model  simulated  plant  growth  and  succession 
on  an  estuarine  salt  marsh,  based  on  an  analysis  of  field  studies  of 


Spartina  alt erni flora,  £3.  patens ,  and  Distichlis  spicata  and  their  re¬ 
sponse  to  tidal  inundations.  The  model  showed  that  plant  growth  and 
succession  across  a  transect  were  controlled  primarily  by  tidal  inunda¬ 
tion,  with  influences  felt  from  radiation  and  temperature,  but  little 
influence  due  to  salinity.  The  authors  state  that  the  model  is  limited 
in  application  due  to  incomplete  parameterization  of  some  of  the  pro¬ 
cesses,  but  has  potential  for  wide  applicability  once  these  processes 
are  better  understood. 

73.  Extensive  mangrove  swamps  are  found  in  tropical  and  subtropi¬ 
cal  systems.  Their  greatest  development  in  this  country  has  occurred  in 
Florida.  These  mangrove  systems  serve  a  similar  function  and  importance 
as  the  salt  marsh  systems  found  in  other  parts  of  the  country.  Lugo, 
Sell,  and  Snedaker  (1973)  first  constructed  a  conceptual  model  of  these 
Florida  mangrove  systems  to  be  used  as  a  tool  in  designing  research  in 
the  mangrove  system.  Inputs  from  this  research  were  later  used  to  pro¬ 
vide  a  simulation  of  the  system. 

71*.  The  simulation  model  was  much  less  complex  than  the  conceptual 
model  and  dealt  only  with  mangrove  production.  As  such,  only  three  com¬ 
partments  were  included  in  the  simulation;  these  were  mangrove  biomass, 
detritus,  and  nutrients.  Rate  coefficients  used  in  the  simulation  were 
based  on  average  and/or  maximum  measured  values.  In  addition,  inputs  to 
the  nutrient  pool  from  other  than  recycling  were  lumped  to  produce  a 
single  number.  The  model,  then,  is  sufficient  for  simulating  the  system 
over  the  short  term  (months,  weeks,  or  days)  without  redefining  the 
constants.  Therefore,  the  usefulness  of  the  model  would  be  a  function 
of  the  long-  or  short-term  nature  of  the  questions  the  model  needed 
to  answer. 

75*  Several  things  are  common  in  all  of  the  models  discussed 
above  and  should  be  mentioned.  First  is  the  multidisciplinary  nature  of 
the  models:  all  incorporate  various  degrees  of  physical,  chemical,  and 
biological  parameterization.  Although  in  many  cases  it  may  be  possible 
to  model  physical  and  chemical  processes  without  significant  biological 
inputs,  it  is  not  possible  to  model  biological  systems  without  these 
other  processes.  Tlerefore,  any  ecosystem  modeling  efforts  should 
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involve  the  efforts  of  experts  in  the  physical,  chemical,  and  biological 
sciences . 

76.  Second,  initial  model  building  for  most  of  the  models  in¬ 
cluded  the  construction  of  a  conceptual  model.  Third,  based  upon  this 
conceptual  model,  a  research  program  was  designed  to  satisfy  the  data 
needs  identified  in  the  conceptual  model.  Lastly,  based  on  the  data 
collected,  the  conceptual  model  was  generally  simplified  and  developed 
into  a  mathematical  simulation  model.  The  need  to  optimize  research  is 
always  important  and,  because  of  the  complexity  of  marshes,  it  is  gener¬ 
ally  not  feasible  to  model  or  measure  all  components  of  these  ecosystems. 
This  is  the  reason  for  the  generalization  and  aggregation  of  ecosystem 
components  seen  in  previous  modeling  efforts.  Any  modeling  efforts  in 
the  future  should,  therefore,  have  well-defined  objectives  sc  that  the 
model  can  provide  answers  for  those  questions  for  which  it  is  designed. 


Relationship  between  the  Marsh  and  the  Estuar\ 


77.  Perhaps  the  most  important  question  facing  the  estuarine 
ecologist  is  the  relation  that  exists  between  the  marsh  and  the  estuary. 
Gunnison  (1978)  has  attempted  to  show  the  relationships  that  exist  in 
his  conceptual  models  of  the  marsh-estuarine  system.  Certainly,  it  can 
be  said  that  the  marsh-estuarine  complex  is  intimately  linked  to  produce 
a  biologically  productive  system.  It  is  generally  accepted  that  the 
major  importance  of  the  estuary  with  respect  to  the  marsh  includes  tidal 
inundation  and  nutrient  input  (both  dissolved  and  particulate);  the 
marsh  is  a  major  contributor  of  detritus  and,  to  a  lesser  extent,  dis¬ 
solved  minerals  to  the  estuary  (Gunnison  1978).  Recently,  however,  this 
concept  of  a  net  import  of  materials  on  or  off  the  marsh  has  begun  to  be 
questioned.  In  the  studies  of  Nixon  et  al.  (1976)  on  a  Rhode  Island 
salt  marsh,  the  inorganic  nutrients  appeared  to  cycle  within  the 
sediment-detritus  system  of  the  embayment  bottom  suggesting  little  export 
of  nutrients  from  the  marsh  during  tidal  exchanges  with  the  estuary. 

78.  The  investigation  of  the  advective  flux  of  material  in  shal¬ 
low  tidal  creeks  has  proved  to  be  very  difficult.  The  usual  method  is 
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to  use  current  meters  and  tide  staffs  and  sample  for  the  material  con¬ 
centration.  Boon  (1978)  has  given  an  illuminating  analysis  of  the  er¬ 
rors  involved  in  calculating  the  flux  of  material  in  a  Virginia  tidal 
creek  using  a  dense  array  of  current  meters  that  record  current  speed 
every  minute.  He  showed  that  substantial  errors  in  calculating  a  net 
flux  are  likely  if  only  a  few  current  meters  are  used  with  sample  times 
of  30  min.  Boon  only  sampled  over  a  number  of  single  tidal  periods. 
Recently,  Baylis-Smith  et  al.  (1979)  in  a  study  of  tidal  creek  flows  in 
an  English  salt  marsh  showed  that  the  net  flux  per  tidal  period  can  vary 
considerably  over  the  spring-neap  cycle;  therefore,  much  longer  sampling 
periods  may  be  necessary  in  the  future.  Finally,  Imberger  (1979)  has 
taken  a  new  approach  to  the  investigation  of  tidal  creek  flows  by  follow¬ 
ing  water  masses  rather  than  measuring  currents.  He  shows  that  the 
water  masses  maintain  their  integrity  over  many  tidal  cycles,  implying 
that  the  exchange  of  detritus  between  the  marsh  and  the  tidal  creek  and 
the  tidal  creek  and  the  estuary  is  probably  small.  There  is  the  possi¬ 
bility,  yet  to  be  investigated,  that  storm  events,  such  as  thunderstorms 
and  storm  surges,  effect  a  large  flux  of  detritus  from  the  marsh  to  the 
estuary  at  infrequent  intervals. 

79-  The  above  discussion  serves  to  point  out  the  difficulty 
alluded  to  in  earlier  discussions  on  modeling  of  estuarine  and  marsh  eco¬ 
systems,  i.e.,  a  lack  of  a  basic  understanding  of  processes  within  these 
systems.  Gardner  and  Kitchens  (1978)  have  reviewed  the  literature  deal¬ 
ing  with  the  exchange  of  sediments  and  chemicals  between  salt  marshes 
and  coastal  waters.  Their  review  showed  that  no  clear  patcerns  of  marsh 
behavior  are  evident  with  the  existing  data. 

80.  Attempts  to  understand  the  marsh-estuarine  complex  are  com¬ 
plicated  by  the  variation  that  exists  within  and  between  even  similar 
marsh  systems  (Gardner  and  Kitchens  1978).  In  a  comparative  evaluation 
of  ten  intertidal  salt  marshes  in  Rhode  Island,  Oviatt,  Nixon,  and 
Garber  (1977)  found  an  apparent  lack  of  intercorrelation  among  the  param¬ 
eters  which  they  studied  (e.g.  marshes  with  large  crab  populations  may 
have  few  fish,  those  with  many  fish  may  not  have  many  shrimp,  etc.). 
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Similarly,  Spartina  production  was  not  always  directly  related  to  animal 
numbers . 

8l.  As  the  earlier  discussions  have  shown,  attempts  to  model  the 
marsh-estuarine  ecosystem  have  generally  emphasized  a  particular  compo¬ 
nent  of  this  system  (i.e.,  either  the  marshes  or  the  estuary).  The 
present  lack  of  understanding  of  the  processes  working  within  this  system 
(i.e.,  detrital  export  and  import),  and  the  relationships  that  exist 
there  (e.g.,  macrophyte  production  and  animal  productivity),  presently 
limit  the  ability  to  model  the  total  system.  Gardner  and  Kitchens  (1978) 
suggest  long-term,  continuous  collection  of  transport  data  in  conjunc¬ 
tion  with  intensive  studies  of  the  intertidal  processes  and  external 
conditions  that  control  exchanges  between  the  marsh  and  the  estuary  to 
alleviate  this  lack  of  knowledge. 
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Representative 


Table  1 


Types  of  Hydrodynamic  Models  of  Estuaries 


Reference 

Dimensions 

Solution  Technique 

Variables 

Equation? 

Forcing 

Elliott  !IKi’5) 

x 

(2  layer) 

Noncoupled  kinematic 
and  dynamic  box  model 

au,l(x) 

uu,V> 

Continuity 

Salt  balance 

River flow 

VS, 

u-momentum 
( approximated ) 

Hamilton  i  975) 
Blumberg  ( *  -75 ) 

Wang  and  Kr  /iz  11979) 

x ,  z ,  t 

Explicit 

Finite  differences 
Semi-implicit 

C(x,t) 

u(x,z,t) 

w(x,z,t) 

s(x,z,t) 

Continuity 

u-momentum 

Salt  conservation 

River flow 
and  tide 
+  wind 

Hodgins  (1979) 

x,t 

(2  layer) 

Finite  differences 
Semi-implicit 

t(x,t ) 
h(x,t ) 

vV*41 

su,l<x,t> 

Continuity 

u-momentum 

Salt  conservation 

Wind 

River flow 

Johns  (1978) 

x  ,z  ,t 

Finite  differences 
Explicit  x-coord 
Semi-implicit 
z-coord 

C(x,t) 

u(x,z,t) 

w(x,z,t) 

E(x,z,t) 

Continuity 
u-momentum 
Turbulent  energy 

Tide 

Reid  et  al.  ( 1977a) 

x,y,t 

Explicit 

Finite  differences 
(subgrid  scale) 
(channels  &  barriers) 

C(x.y.t) 

tr(x,y,t) 

^(x.y.t) 

Continuity 

TT-momentum 

V-moraentum 

Tide 
&  wind 

Pearson  and 

Winter  (1977) 

x,y,t 

Finite  elements 

(time  decomposed  into 
tidal  harmonies  w) 

t(x,y.“) 

u(x,y,u) 

V(x,y,u) 

Continuity 

u-momentum 

V-momentum 

Tide 

Liu  and 

Leendertse  (1978) 

x.y.z.t 

Finite  differences 
Semi-implicit 

C(x,y,t) 
u(x,y,z,t ) 
v(x,y,z,t ) 
w(x,y,z,t) 
s(x,y,z,t ) 
E(x,y,z,t ) 

Continuity 

u-momentum 

v-momentum 

Salt  conservation 
Turbulent  energy 

Tide 

Wind 

Riverflow 

Variable 

x 

y 

2 

t 


S 

E 

f. 

h 

VS, 

Superscripts  utl 


_ Definition _ 

Longitudinal  horizontal  axis 
Across-channel  horizontal  axis 
Depth  axis 
Time 

Harmonic  decomposition  of  time  (i.e*  <(x,t) 

*  EA  (x)exp(iw  t) ) 
n  n 

Velocities  in  the  x-,  y-,  and  z-directions,  respectively 

Depth -ayeraged  velocities  in  the  x-  and  y-directions, 
respectively 

Salinity 

Turbulent  kinetic  energy 

Free  surface  elevation 

Interface  depth  of  two-layer  models 

Vertical  and  horizontal  turbulent  diffusion  coeffi¬ 
cients,  respectively 

Upper  and  lower  layers,  respectively 


Application 
Potomac  Estuary 


Rotterdam  Waterway 
Potomac  Estuary 


Alberni  Inlet , 
British  Columbia 


Rectangular 

uniform  channel 


Sabine-Calcasieu 
area  of  gulf 
coast 


Hood  Canal  Wash. 
(Jamart  personel 
communication) 

Long  Island  Sound 
San  Francisco  Bay 

Chesapeake  Bay 


Table  2 


Reference 


Barrett  and 
Mollowney  (1972) 


Nihoul  et  al. 
(1979) 


Leendertse  and 
Gritton  (1971) 
Leendertse  and 
Liu  (1975) 

Chen  and  Orlob 
(1972) 


DiToro  et  al. 
(1977) 


Dahl-Madsen 

(1978) 


Seip  (1979) 
and  Seip  et  al. 
(1979) 


Representative  Types  cf  Water  Quality  Models  cf  Estuaries 


Hydrodynamic  Circulation 
Submodel* 

Water  Quality 

State  Variables 

Kinetics  and  Comments 

Forcing 

1-D  ti dally  averaged 
box  model 

Organic  carbon,  organic 
nitrogen,  ammonia, 
nitrate,  DO 

Linear 

Anerobic  kinetics  in¬ 
cluded  when  DO  <  5% 
saturation 

Tidal  amplitude 
riverflow,  wind, 
carbonacous  and 
nitrogenous  loads, 
temperature 

1-D  time-dependent 
dynamic  transport 
diapers  ion  model 

DO,  nitrate,  ammonia, 
iron,  manganese 

Thermodynamic  equilibrium 
assumed.  Redox  poten¬ 
tials  control  reaction  ’ 
rates  assumed  to  be  due 
tc  bacterial  activity 
Highly  polluted  -  no 
nutrient  limitation 

Tide,  riverflow 

2-D  (horizontal) 
t ime-dependent 
dynamic  transport- 
dispersion  model 

Conforms,  salinity  (chlo¬ 
rides),  DO-BOD 

Linear 

Tide,  fresh  and 
storm  drain  flow, 
BOD  load,  wind 

1-D  time-dependent 
link-mode  dynamic 
transport-di spersi on 
model 

Temperature,  salinity 
ammonia  nitrogen, 
nitrate  nitrogen, 
phosphorus,  DO-BOD, 
algae  (2),  zooplankton, 
fish  (2) 

Linear  and  Monod 

Tide,  riverflow, 
waste  and  nutrient 
loads,  wind,  light 

1-D  tldally  averaged 
transport-dispersion 
model 

DO-BOD,  phosphorus  (or¬ 
ganic  and  inorganic), 
silica,  nitrogen  (ni¬ 
trate,  ammonia,  organic), 
phytoplankton,  zooplank¬ 
ton 

Linear  and  Monod 

Tide,  riverflow, 
wind,  temperature, 
light ,  carbonacous 
and  nitrogenous 
waste  loads 

Box  models  or  1-D  time- 
dependent  transpert- 
dispersion  model 

Phytoplankton  (carbon,  ni¬ 
trogen,  phosphorus),  zoo¬ 
plankton  (carbon),  detri¬ 
tus  (carbon,  nitrogen, 
phosphorus),  inorganic 
nitrogen,  inorganic  phos- 
phosphorus,  DO  -  BOD, 
sediment  carbon 

Linear  and  Monod 

Riverflow,  tide, 
temperature, 
light,  nutrient 
and  organic  loads 

None  (observed  hydro- 
graphic  regime  used 
as  forcing  for  four- 
layer,  four-  horizon¬ 
tal  section  model) 

Benthic  algae  (8  age 

classes),  biomass,  zinc, 
iron 

Model  of  population 

dynamics  with  toxicity- 
related  mortality 

Light,  salinity, 
temperature,  ni¬ 
trogen  (limiting 
nutrient ) 

_ Application 

Thames  Estuary 
(U.K. } 


Scheldt 

(The  Netherlands ) 


Jamaica  Bay 
(N.Y.) 


San  Franc i see  Bay 
Delta  System 


San  Franc Ibcc 
Bay-Delta  System 
Fctcraac  Estuary 


Various  Danish  fjords 
and  estuaries 
Principally 
Limfjorden 


Hardanger  fjord 
Sorfjorden 
Trondheims  fjorden 
( Norway ) 


dimensional j  2-D  •  tec  dimensional. 
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